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The redshift evolution of ACDM halo parameters: 
concentration, spin and shape 
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ABSTRACT 

We present a detailed study of the redshift evolution of dark matter halo structural 
parameters in a ACDM cosmology. We study the mass and redshift dependence of 
the concentration, shape and spin parameter in N-body simulations spanning masses 
from 10^° h~^ Mq to 10^^ h~^ Mq and redshifts from to 2. We present a series of 
fitting formulae that accurately describe the time evolution of the concentration-mass 
(cvir — Mvir) relation since z = 2. Using arguments based on the spherical collapse 
model we study the behaviour of the scale length of the density profile during the as- 
sembly history of haloes, obtaining physical insights on the origin of the observed time 
evolution of the Cvir — Mvir relation. We also investigate the evolution with redshift of 
dark matter halo shape and its dependence on mass. Within the studied redshift range 
the relation between halo shape and mass can be well fitted by a redshift dependent 
power law. Finally we show that although for z = the spin parameter is practically 
mass independent, at increasing redshift it shows a increasing correlation with mass. 



Key words: galaxies: haloes - cosmology:thcory, dark matter, gravitation - 
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1 INTRODUCTION 

Observational evidence (e.g. Spergel et al. 2007; Komatsu 
et al. 2009) favours the hierarchical growth of structures 
in a universe dominated by cold dark matter (CDM) and 
dark energy (A), the so called ACDM universe. Within this 
paradigm dark matter collapses first into small haloes which 
accrete matter and merge to form progressively larger haloes 
over time. Galaxies are thought to form out of gas which 
cools and collapses to the centres of these dark matter haloes 
(e.g. White & Rees 1978). In this picture the properties of 
galaxies are expected to be strongly related to the properties 
of the dark matter haloes in which they are embedded (e.g. 
Mo, Mao & White 1998, Button et al. 2007). 

It has been shown by several studies that the struc- 
tural properties of dark matter (DM) haloes are dependent 
on halo mass: for example higher mass haloes are less con- 
centrated (Navarro, Frank & White 1997, hereafter NFW; 
Eke et al. 2001; Bullock et al. 2001a; Kuhlen et al. 2005; 
Maccio et al. 2007; Neto et al. 2007; Gao et al. 2008; Duffy 
et al. 2008; Maccio, Dutton, & van den Bosch 2008 (hereafter 
M08); Klypin et al. 2010), and are more prolate (Jing & Suto 



* Email : jcmunoz@aip.de 
t CITA National Fellow 



2002; AUgood et al. 2006; Gottlober & Yepes 2007; Bett et 
al. 2007; Maccio et al. 2007; M08) on average. The situation 
is less clear for the spin parameter, at z = there seems to 
be no mass dependence (Maccio et al. 2007; M08) or at least 
a very weak one (Bett et al. 2007), while for increasing val- 
ues of the redshift a possible mild correlation between spin 
and mass seems to arise (Knebe & Power 2008). 

In M08 the properties of DM haloes were studied in 
ACDM universes whose parameters were fixed by the one, 
three and five-year release of the WMAP mission (WMAP5; 
Komatsu et al. 2009). In that study the attention has been 
focused on the structural parameters of virialized haloes 
and their correlations at the present epoch, z — 0. In this 
work we extend this previous analysis to higher redshifts and 
study how the scaling relations of DM haloes change with 
time. 

As in M08 we use a large suite of N-body simula- 
tions in a WMAP5 cosmology with different box sizes to 
cover the entire halo mass range relevant for galaxy forma- 
tion: from lO^^/i"^ Mq (haloes that host dwarf galaxies) to 
10^^h~^ Mq (massive clusters). We use these simulations to 
investigate the evolution of concentrations, spin parameter 
and shapes of dark matter haloes through cosmic time. 

Similar studies have been already conducted in the past, 
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mainly using lower numerical resolution and/or a smaller 
mass range (but with few recent exceptions). 

Navarro, Frank & White (1997) proposed that the char- 
acteristic density of dark matter haloes was directly propor- 
tional to the density of the universe at time of formation, 
making possible to connect today properties of the dark mat- 
ter density profile to the halo formation history and to the 
evolution of the expanding universe. This idea was then ex- 
panded by Wechsler et al. (2002), who found a clear con- 
nection between the mass growth of dark matter haloes and 
the definition of the formation time; connecting directly the 
growth history of dark matter haloes to the evolution of 
their concentration parameter. 

In a series of papers, Zhao et al. (2003a, 2003b, 2009) 
have re-addressed the problem of the evolution of dark 
matter halo density profile and the mass accretion history. 
Zhao's main result was the finding of a correlation between 
Ts and the characteristic mass of dark matter haloes, Ms, de- 
fined as the mass inside r^. Thanks to this correlation they 
were able to model the time evolution of the concentration 
parameter in a cosmology- free fashion. 

A comparison of the difi'erent approaches to predict 
dark matter halo concentrations was performed by Neto et 
al. (2007). They made a detailed comparison of the Wech- 
sler et al. (2002) and Zhao et al. (2003) models for the time 
evolution of dark matter halo masses and their resulting 
predictions for halo concentration. Neto et al. (2007) found 
that although these models could match the average concen- 
tration reasonably well, they performed very poorly in many 
cases, because their models for halo mass evolution were not 
able to satisfactorily reproduce "real" mass growth histories 
from N-body simulations. 

However, the evolution of dark matter halo properties 
does not reduce to the concentration parameter only; halo 
shape and spin parameter are also important quantities that 
could influence the properties of the hosted galaxy. All- 
good et al. (2006) studied the mass, radius, redshift and 
cosmology dependence (via variations of ag,) of the shape 
of dark matter haloes; while the environment dependence 
of the shape has been addressed in Hahn et al. (2007a,b). 
The distribution of the spin parameter of dark matter 
haloes has been studied in several works (e.g. Bullock et 
al. 2001b; Maccio et al. 2007; Bett et al. 2007; M08; Knebe 
& Power 2008; Davis & Natarajan 2009; Antonuccio-Delogu 
et al. 2010) as well as the correlation between halo angular 
momentum and large scale structure (e.g. Bailin & Stein- 
metz 2005, Bett et al. 2010). One of the conclusions was 
that the correlation between spin and mass (almost absent 
at 2 = 0) seems to increase with increasing redshift. Such a 
behavior could have important influences in the modeling of 
properties of galaxies at high redshift. 

Although it has been shown that the inclusion of bary- 
onic physics may affect the properties of the dark matter 
distribution (e.g. Gnedin et al. 2004, Kazantzidis et al. 2004, 
Libeskind et al. 2010, Knebe et al. 2010, Bett et al. 2010), it 
also known that the strength of this effect strongly depends 
on the implemented baryonic physics (Duffy et al. 2010). 
The aim of this work is to extend the findings of M08 to 
higher redshift and to provide a comprehensive study of the 
evolution and correlation of the properties of DM haloes 
from z = 2 to the present day. Such a study is a basic ingre- 
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Table 1. Table of simulations used in this work. Note that the 
name of the simulation is related to the box size in units of Mpc. 
N represents the number of total particles in the box. e represents 
the force softening length in units of kpc h~^ and the last column 
gives the number of haloes with more than 500 particles at z = 
and 2 = 2. Masses of particles are in units of h~^ Mq and box 
sizes in units of Mpc h^^, with h = 0.72. 



dient to understand the complex problem of galaxy forma- 
tion and evolution. 

Our paper is organized as follows. In Section [2] we de- 
scribe the set of simulations we used, present definitions and 
describe the numerical procedures. In Section [3] we present 
our results for the mass and redshift dependence of the con- 
centration parameter, present a novel and physically mo- 
tivated approach to understand the time evolution of the 
concentration parameter following arguments based on the 
spherical collapse model. Next we discuss our results for the 
mass and redshift dependence of the shape (Section UJ and 
spin parameter (Section [Sjl. Finally in Section |6] we summa- 
rize our findings and discuss on our results and their influ- 
ence on the modeling of galaxy properties. 



2 NUMERICAL SIMULATIONS 

All simulations have been performed with pkdgrav, a tree 
code written by Joachim Stadel and Thomas Quinn (Stadel 
2001). The code uses spline kernel softening, for which the 
forces become completely Newtonian at 2 softening lengths. 
Individual time steps for each particle are chosen propor- 
tional to the square root of the softening length, e, over the 
acceleration, a: Ati = r]^ tjai. Throughout, we set r\ — 0.2, 
and we keep the value of the softening length constant in 
comoving coordinates during each run. The physical values 
of e at 2 = are listed in Table[T] Forces are computed using 
terms up to hexadecapole order and a node-opening angle 
Q which we change from 0.55 initially to 0.7 at 2: = 2. This 
allows a higher force accuracy when the mass distribution 
is nearly smooth and the relative force errors can be large. 
The initial conditions are generated with the GRAFIC2 pack- 
age (Bertschinger 2001). The starting redshifts Zi are set to 
the time when the standard deviation of the smallest den- 
sity fluctuations resolved within the simulation box reaches 
0.2 (the smallest scale resolved within the initial conditions 
is defined as twice the intra-particle distance). 

Table [T] lists all of the simulations used in this work. 
We have run simulations for several different box sizes, which 
allows us to probe halo masses covering the entire range 10^" 
h~^ Mq < M < 10^^ h~^ Mq. In addition, in some cases we 
have run multiple simulations for the same cosmology and 
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box size, in order to test for the impact of cosmic variance 
(and to increase the final number of dark matter haloes). 

In all of the simulations, dark matter haloes are iden- 
tified using a spherical overdensity (SO) algorithm. We use 
a time varying virial density contrast determined using the 
fitting formula presented in Bryan and Norman (1998). We 
include in the halo catalogue all the haloes with more than 
500 particles inside the virial radius (A'vir > 500). 

We have set the cosmological parameters according to 
the fifth-year results of the Wilkinson Microwave Anisotropy 
Probe mission (WMAP5, Komatsu et al. 2009), namely, 
Qm = 0.258, Ql = 0.742, n = 0.963, h = 0.72, and 
(78 = 0.796, where ilm and Ol are the values of the den- 
sity parameters at z=0. 



2.1 Halo parameters 

For each SO halo in our sample we determine a set of param- 
eters, including the virial mass and radius, the concentration 
parameter, the angular momentum, the spin parameter and 
axis ratios (shape). Below we briefly describe how these pa- 
rameters are defined and determined. A more detailed dis- 
cussion can be found in Maccio et al. (2007, 2008). 



2.1.1 Concentration parameter 

To compute the concentration of a halo we first determine 
its density profile. The halo centre is defined as the location 
of the most bound halo particle (we define the most bound 
particle as the particle with the lowest potential energy, no 
care about binding energy is taken here), and we compute 
the density (pi) in 50 spherical shells, spaced equally in log- 
arithmic radius. Errors on the density are computed from 
the Poisson noise due to the finite number of particles in 
each mass shell. The resulting density profile is fitted with 
a NFW profile: 



PJr) ^ ^c 

pc (?-/r-s)(H-r/rs)2' 



(1) 



During the fitting procedure we treat both r^ and 5c as 
free parameters. Their values, and associated uncertainties, 
are obtained via a x^ minimization procedure using the Lev- 
enberg & Marquardt method. We define the r.m.s. of the fit 



N 
Orms = -r7 7 (Inp, - Inpni)^ 



N ^^ 



(2) 



where pm is the fitted NFW density distributions Finally, 
we define the concentration of the halo, Cvir = Rvii/rs, using 
the virial radius obtained from the SO algorithm, and we 
define the error on logc as {ar^/rs)/ln{10), where ar^ is the 
fitting uncertainty on r^. 



^ A more conservative notation for the r.m.s of the fit would be 
cr'i. Nevertheless we keep the use of prms to be consistent with 
the notation used in Maccio et al. (2007, 2008) 



2.1.2 Spin parameter 

The spin parameter is a dimensionless measure of the 
amount of rotation of a dark matter halo. We use the defi- 
nition introduced by Bullock et al. (2001b): 

Jvir 



A 



(3) 



V 2Aivir VvirJ^vir 

where Mvir is the mass interior to -Rvir, Jvh is the total 
angular momentum of that mass distribution and Vvir is its 
circular velocity at the virial radius. See Maccio et al. (2007) 
for a detailed discussion and for a comparison of the different 
definitions of the spin parameter. 

2.1.3 Shape parameter 

Determining the shape of a three-dimensional distribution 
of particles is a non-trivial task (e.g., Jing & Suto 2002). 
Following AUgood et al. (2006), we determine the shapes of 
our haloes starting from the inertia tensor. As a first step, we 
compute the halo's 3x3 inertia tensor using all the particles 
within the virial radius. Next, we diagonalize the inertia 
tensor and rotate the particle distribution according to the 
eigenvectors. In this new frame (in which the moment of 
inertia tensor is diagonal) the ratios as/ai and 03/02 (where 
di ^ 02 ^ 03) are given by: 



03 
ai 



'E"^i2f 



E' 



03 

a2 



'E' 



E '^iVf 



(4) 



Next we again compute the inertia tensor, but this time 
only using the particles inside the ellipsoid defined by ai, a2, 
and 03. When deforming the ellipsoidal volume of the halo, 
we keep the longest axis (ai) equal to the original radius of 
the spherical volume (i?vir). We iterate this procedure until 
we converge to a stable set of axis ratios. Although this iter- 
ative procedure can indeed change the mass contained inside 
of the ellipsoid, we checked that variations are nevertheless 
below 20%. We will therefore not consider those changes in 
mass when showing mass-shape relations, and we will work 
always with virial masses. 



2.2 Relaxed - Unrelaxed haloes 

Our halo finder (and halo finders in general) does not distin- 
guish between relaxed and unrelaxed haloes. There are many 
reasons why we might want to remove unrelaxed haloes. 
First and foremost, unrelaxed haloes often have poorly de- 
fined centers, which makes the determination of a radial den- 
sity profile, and hence of the concentration parameter, an 
ill-defined problem. Moreover unrelaxed haloes often have 
shapes that are not adequately described by an ellipsoid, 
making our shape parameters ill-defined as well. 

Following Maccio et al. (2007), we decide to use a combi- 
nation of two different parameters prms and Xoff to determine 
the dynamical status of a given dark matter halo. The first 
quantity prms is defined as the r.m.s. of the NFW fit to the 
density profile (performed to compute Cvir). While it is true 
that Prms is typically high for unrelaxed haloes, haloes with 
relatively few particles also have a high prms (due to Pois- 
son noise) even when they are relaxed; furthermore, since 
the spherical averaging used to compute the density profiles 
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has a smoothing effect, not ail unreiaxed haioes have a high 
Prms- In order to circumvent these probiems, we combine the 
value of prms with the XoS parameter, defined as the distance 
between the most bound particle (used as the center for the 
density profile) and the center of mass of the halo, in units 
of the virial radius. This offset is a measure for the extent 
to which the halo is relaxed: relaxed haloes in equilibrium 
will have a smooth, radially symmetric density distribution, 
and thus an offset that is virtually equal to zero. Unreiaxed 
haloes, such as those that have only recently experienced 
a major merger, are likely to reveal a strongly asymmetric 
mass distribution, and thus a relatively large XoH- Although 
some unreiaxed haloes may have a small XoS, the advan- 
tage of this parameter over, for example, the actual virial 
ratio, 2T/V , as a function of radius (e.g. Maccio, Murante 
& Bonometto 2003), is that the former is trivial to evaluate. 
Following Maccio et al. (2007), we split our halo sample into 
unreiaxed and relaxed haloes. The latter are defined as the 
haloes with prms < 0.5 and Xo« < 0.07. About 70% of the 
haloes in our sample qualify as relaxed haloes a.t z = Q. 

To check for the effect of changing the definition of re- 
laxed haloes, we have computed the median concentration 
(as shown in the next section) using different values of Xofi. 
Changing the value of this parameter by 25% (above and 
below 0.07) induces changes no larger than 5% in the me- 
dian concentration of dark matter haloes. We conclude that 
choosing Soff = 0.07 our results are robust enough against 
variations in the definition of relaxed population of haloes. 
In what follows we will just present results for haloes which 
qualify as relaxed. 



CONCENTRATION: 
DEPENDENCE 



MASS AND REDSHIFT 



In figure[T]we show the median Cvir— Mvir relation for relaxed 
haloes in our sample at different redshifts. Haloes have been 
binned in mass bins of 0.4 dex width, the median concentra- 
tion in each bin has been computed taking into account the 
error associated to the concentration value (see 12.1.11 and 
M08). In our mass range the Cvir — Mvir relation is well fitted 
by a single power law at almost all redshifts. Only for 2 = 2 
we see an indication that the linearity of the relation in log 
space seems to break, in agreement with recent findings by 
Klypin et al. (2010). 

The best fitting power law can be written as: 

log(c) = a{z) log(Afvir/[/i"' Mo]) + b(z) (5) 

The fitting parameters a{z) and b[z) are functions of 
redshift as shown in figure (2] The evolution of a and 6 can 
be itself fitted with two simple formulas that allow to recon- 
struct the Cvir — Mvir relation at any redshifts: 



a{z) = wz 
biz) 



+ 



P 



(6) 

(7) 



(z + 'y) ' (z + j)^ 

Where the additional (constant) fitting parameters have 
been set equal to: w = 0.029, m = 0.097, a = -110.001, 
l3 = 2469.720 and 7 = 16.885. Figure [3] shows the recon- 
struction of the Cvir — Mvir relation for different mass bins as 
a function of redshift using the approach described above. It 
shows that our (double) fitting formulas are able to recover 
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Figure 1. Mass and redshift dependence of the concentration 
parameter. The points show the median of the concentration as 
computed from the simulations, averaged for each mass bin. Lines 
show their respective linear fitting to eq. [5] 



the original values of the halo concentration with a precision 
of 5%, for the whole range of masses and redshifts inspected. 
It has been shown by Trenti et al. (2010) that using A'vir be- 
tween 100 and 400 particles is enough to get good estimates 
for the properties of halos, nevertheless in order to look for 
systematics we re-computed Cvir varying the minimum num- 
ber of particle inside -Rvir, using 200, 500 and 1000 particles. 
No appreciable differences (less than 2%) were found in our 
results for the median. We also checked that our results do 
not change notably by changing the definition of "relaxed" 
haloes (i.e. changing the cut in XoS or prms). 

It is interesting to compare our results with M08, which 
shares some of the simulations presented in this work. Our 
results for the Cvir — Mvir relation at 2; = are slightly dif- 
ferent to those presented in M08: we found a = —0.097 and 
b = 2.155, while M08 found a = -0.094 and b = 2.099. 
The difference is less than 3%, and it is mainly due to low 
mass haloes. In this work we included three new simulations, 
B30, B90 and B3OO2. Two of these (B30 and B90) increase 
the statistics of our halo catalogs at the low mass end, pro- 
viding a better determination of the Cvir — Mvir relation for 
M « 10^^/i~^ Mq. We are confident that the inclusion those 
new simulations led to an improvement over the results of 
previous works. 



3.1 Understanding the concentration evolution 

We want now to explore in more detail the physical mecha- 
nism driving the mass and redshift dependence of the con- 
centration parameter, this understanding could take us to a 
better interpretation of the time evolution of the dark mat- 
ter density profile. As Cvir is defined as the ratio between R^ir 
and Ts we will look at the time evolution of those quantities 
for different halo masses, and try to see if we can extract 
some physical insights about the evolution of Cvir. 
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Figure 2. Redshift evolution of the parameters a{z) and b{z) of 
the linear fitting of the Cvir — Mvir relation. Points represent the 
result of the fitting shown in figure [l] for each redshift while lines 
are the fit to the functions in eq: . [6] and [7] Error bars show the 
error obtained from the fitting to eq. \5\ 
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Figure 3. Reconstruction of the Cvir — Mvir relation for different 
mass bins as a function of redshift as computed from the model 
fitting function, equations [5] to [7] Points show the data of our 
simulations while the lines show the result of the reconstruction. 
Error bars show the error in the median concentration for each 
redshift. 



Since the evolution of the concentration parameter is 
strongly correlated to the mass growth history of the dark 
matter halo (e.g. Wechsler et al. 2003; Zhao et al. 2009) we 
need to construct merger trees for haloes in our simulation 
boxes. Here we briefly describe the procedure used to build 



merger trees (for more details see Neistein, Maccio & Dekel 
2010). 

We link haloes A and B, with particle number Na and 
Nb, in two consecutive snapshots at redshift za and zb 
{za < Zb), if they success to satisfy a list of requirements: 

• if Nb < Na, they have at least 0.57Vs particles in com- 
mon 

• if Na < Nb, they have at least 0.5A''a particles in com- 
mon 

• Halo B does not contribute to any other halo in za with 
more particles than it does for halo A 

We will assume that the evolution of the structural 
properties of the halo is traced by the most massive pro- 
genitor, therefore we will present results only for the main 
branch of the tree. We follow along the tree the evolution of 
Afvir, -Rvir, ^g aud Cvir- Siucc individual histories of haloes 
can be very different we focus our study on the time evolu- 
tion of the mean value of these quantities. These averages 
are computed binning the histories by mass at 2 = O with 
each mass bin having width of 0.4 dex. One may argue that 
following only the merger tree main branch of haloes inden- 
tified at z=0 could introduce a bias when comparing to the 
evolution of haloes at earlier cosmic times. In order to verify 
the stability of our results we rebuilt the merger histories 
for all haloes starting from z=0.5 and z=1.0 (i.e. without 
taking into account the future evolution of an halo, like, for 
example, being or not in the final catalogue at z=0). We 
found no differences in the evolution of the mean values for 
-Rvir or Afvir with respect to the original merger tree built 
from z=0. 

For the averages we include only histories with the same 
length, that is, histories that in the entire box evolve from 
the same Zinit to z — 0. Zi„it was chosen by numerical- 
statistic reasons. As we are computing the mean on several 
properties of haloes in the merger histories along the time, 
we wanted to make sure to have at each redshift a statisti- 
cally meaningful number of histories. The extension in red- 
shift is then constrained by the resolution of our simulation 
and the number of histories we choose for our study. At the 
end we found Zinit ~ 3.5 to be a good compromise. More- 
over we only use histories in which the fitting of the density 
profile was successful in the majority of the snapshots. Any 
halo was allowed to be "unrelaxed" for a maximum of 2 con- 
secutive snapshots. In this case, we used as criterion for the 
goodness of the density profile fit prms > 0.5. Using both 
our criteria of relaxation (prms and a::oft) the number of his- 
tories becomes too small to be statistically meaningful. To 
be able to follow the histories longer in redshift we reduced 
the minimum number of particles per halo (inside the virial 
radius) from 500 (as used in the previous section) to 200. 
We verified that using 200 particles inside the virial radius 
the fitting of the density profile is still acceptable for the 
purposes of this section. 

Finally, since we are constraining the mass histories 
in redshift extent, we cannot overlap data among different 
boxes, then we have computed average quantities for indi- 
vidual boxes independently. In the following we will present 
results coming from the B90 box, other boxes show very 
similar behaviours. Applying all of this selection criteria on 
B90 results in a set of 2300 histories with the different mass 
bins having between 35 and 800 halos each. 
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Figure 4. Mean mass accretion history for haloes in the 90 Mpc 
box computed for different mass bins according to their mass at 
z=0 (from 1.68 X 10^^ to 2.8 X lO^^.S ^^-i j^^ ^-^^-^ ^^^^^ ^^^ ^iin 
of width 0.4 dex) as described in the main text. Here is evident the 
different shape of the mass growth history of haloes of different 
masses. Lines show the best fit to Eq[8]which shows to be a good 
description for all of our mass histories. The values of (7, /9) for 
the fit are (0.649,0.273), (0.865,0.259) and (0.900,0.045) for each 
line (low to high mass) respectively. 



In figure |4] we show the average mass accretion history 
(normalized to the z — value) for dark matter haloes in 
three different mass bins centered on: 1.68 x 10^^h~^ Mq, 
1.0 X W^^h'^Me and 6.1 x lO^^/i'^M©. The data are 
well described by the two parameter function (McBride et 
al. 2009) 



M{z) = Mo(l + zf exp(-7z) 



(8) 



Here A'lo is the mass at 2; = and /3 and 7 are free 
parameters related to the mass growth rate at low z. Note 
that although for most haloes /3 7^ 0, when 13 — the pro- 
file assumes the exponential shape adopted by Wechsler et 
al. (2003) with 7 = ln(2)/z/. 

In figure [S] (top panel) we show the redshift evolution 
of the virial radius for haloes with final masses: 1.68 x 
10"/i~^Mq, 1.0 X lO^^/i"^M0 and 6.1 x 10^^/i"^Mq. For 
all mass scales the virial radius grows with decreasing red- 
shift reaching a maximum, and then starts a slow decrease. 
The redshift at which this maximum is reached depends on 
the mass of the halo, lower mass haloes reach that maximum 
earlier than more massive ones. Given the definition of -Rvir: 



RMz) 



3Mvi,(2) 
'inA^ir{z)pc(z) 



1/3 



(9) 



its time evolution can be understood as follows. At high red- 
shift the radius of the virialized region _Rvir grows due to the 
growth of the halo mass M{z) but is also subject to the ef- 
fects of the cosmological background via [Avii{z)pc{z)]~^''^ , 
which in this case is a slowly decreasing function with de- 
creasing redshift. At high redshift the growth of M{z) dom- 
inates, forcing _Rvir to grow, at low redshift the growth rate 



0.4 

0.35 

0.3 



= 1.68x10 'M,- 



= 1.00 X lO^Mrh'' 
=6.10xl0''^M„h'' 



t 


0.25 - 






^ ^ 


0.2 - 


N 










0.15 - 


« 






0.1 - 

r 

0.07 - 




1 — 1 


0.06 - 


j= 


0.05 - 






p. 




^ 


0.04 - 












0.03 - 


;-^ 






0.02 - 




0.01 E 




0.5 



1.5 



2.5 



3.5 



Figure 5. (Top) Time evolution of the virial radius as computed 
from the averaged merger histories for haloes with final mass of 
1.68 X lOli, 1.0 X IO12 and 6.1 x lO^^h-^MQ. (Bottom) Time 
evolution of the averaged scale length Tb of the dark matter haloes 
in the same mass bins. The mass binning is the same as in figure 
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of the halo mass becomes weaker compared to the decrease 
of the factor [Avir(2)pc(2)]~^ , slowing the growth of 7?vir- 
This behaviour depends on the halo mass: at different red- 
shifts haloes of different masses grow at different rates (see 
figure |4ll, on the other hand, for a given redshift, the fac- 
tor [Avir{z)pc{z)]~^'^ is mass independent. As a result the 
point at which 7?vir reaches its maximum happens later in 
time for massive haloes than for low mass ones. 

The bottom panel of figure [5] shows the redshift evolu- 
tion of the halo scale length r^, which has a trend similar to 
-Rvii-: it grows with time, reaches a maximum and then starts 
to decrease. Once again the redshift at which this maximum 
is achieved depends on the mass of the halo, with low mass 
haloes reaching the maximum at earlier times. 

The behaviour of Vs (and i?vir) is strongly reminiscent 
of the behaviour of a perturbation in the spherical col- 
lapse model. It seems to suggest that the inner region of 
an halo (within Ts) evolves in a decoupled way compared 
to the global perturbation (within i?vir). It is then possi- 
ble to model the inner region as a perturbation of density 
ps (defined as the density inside Vs) that evolves within the 
background density Pvir(z) = Avir(z)pc(z)- 

In analogy with the spherical collapse model we want to 
look for the evolution of the density contrast of this pertur- 
bation: As{z) = ps{z)/pvii{z). This inner density contrast is 
well described by the following fornmla: 



Ae(2) 



(10) 



z + e(M) 

where A = 50 and e(M) = 0.3975 log(Mvir (2 = 
0)/[/i~^ M©]) - 4.312. Equation [TOl implies that: i) ps > Pvir 
at all redshifts, ii) Ag is a growing function of the redshift, 
implying a faster growth of the inner density with respect 
to the mean density of the halo and iii) As depends on the 
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final mass of the lialo, and it iias larger values for high mass 
haloes. 

It is now possible to interpret the evolution of the 
Cvir — Mvir relation in the light of our findings. Figure [3] 
shows that the growth rate of the concentration depends on 
the halo mass, with low mass haloes experiencing a faster 
concentration evolution. Moreover, at fixed mass, the evo- 
lution of the Cvir — Mvir relation is faster at lower redshifts. 
As can be seen from figure O at early times r^ grows si- 
multaneously with i?vir, then, with decreasing redshift, the 
growth of Ts slows down, reaches a maximum, and starts to 
decrease. When _Rvir grows together with rg, the concentra- 
tion of the halo stays approximately constant or slightly in- 
creases. Then when Vs slows down its growth rate and starts 
to "contract" the concentration of the halo grows rapidly. In 
the case of high mass haloes the time at which the dynam- 
ics of the inner region decouples from the outermost part, 
happens at later times compared to low mass haloes. As a 
consequence more massive haloes have a more extended pe- 
riod in which the concentration is a slowly growing function 
of time as shown in figure [S] Only in recent cosmic times, 
due to the late collapse of the inner part the concentration 
starts to grow at a higher rate also for larger masses. 

A crucial point in this process is the moment when the 
inner part of the halo decouples from the outer one. This 
specific time is strongly mass dependent with smaller struc- 
tures (i.e. the more non linear ones) decoupling earlier. This 
effect explains the change in the slope of the Cvir — Mvir re- 
lation, as described in figure [5] by the evolution of a{z), 
and also explains why a simple scaling of the redshift zero 
Cvir — Mvir relation with a factor of the form (1 + z)" (e.g. 
Bullock et al. 2001a) is not able to reproduce the simulated 
data. 



4 HALO SHAPES 

We now turn our attention to the evolution of halo shape 
with redshift and its dependence on halo mass. As was 
described in section 12.1.31 for each dark matter halo we 
compute the axis ratios s — a-s/ai and p = a^ja-i^ where 
cii ^ 02 ^ as are the major, intermediate and minor axis 
of the halo mass distribution. It is worth to note that in 
our notation using only s and p is enough to determine the 
shape of a dark matter halo. The condition for a halo to be 
oblate will be s ~ p < 1, while the condition to be prolate 
will be s < p with p ~ 1, where the obvious condition for 
sphericity is s ~ p ~ 1 and triaxiality will be any other not 
filling any of the previous requirements. Figure [S] shows the 
evolution of the shape parameter s as a function of redshift 
and mass. Figure [7] is the analogue of figure [6] but for the 
p parameter, only results for haloes with A'vir > 1000 are 
shown; after several convergence tests we found that this 
number ensures both numerical stability of the axis deter- 
mination and a fairly large statistical sample. Both figures 
show that on average halos are preferentially triaxial, where 
the most massive haloes tend to be the most ellipsoidal ones, 
while lower mass haloes tend to be closer to spherical , and 
this trend seems to be redshift independent. This seems to 
suggest a simple scenario where the most massive haloes are 
the more extended ones (and less concentrated) and hence 
more strongly affected by tidal torques. For each redshift 
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Figure 6. Mass and time evolution of the shape of dark matter 
haloes quantified via the s = a^/ai ratio. Points represent our 
data while the dashed line shows the fitting to equation 1121 to 
^ = 0. Similar results for the fitting are obtained for z > and 
for clarity are not shown in the plot. The color code is the same 
used in figure [T] and shows different redshifts. Values of the fit 
parameters for all redshifts can be found in the appendix. 



we fit the data with a power law of the form of eq: [5] The 
resulting parameters of this fit are listed in the Appendix. 

In order to make a more direct comparison with the 
results of AUgood et al. (2006) we also computed the shape 
parameter s using only particles within the inner 30% of the 
virial radius so. 3 for haloes with A^vir > 4000. The results are 
shown in figure |8] AUgood et al. (2006) proposed a fitting 
formula for the shape parameter given by. 



s{z,A'I) = a 



Mvir 

M.(z) 



(11) 



where M« (z) is the characteristic non-linear mass at z such 
that the rms top-hat smoothed overdensity at scale a{M*, z) 
is 5c ~ 1.68, and a and /3 are free parameters. The quantity 
Mt{z) should in principle contain all the information about 
the cosmological model, making the other fitting parameters 
cosmological independent. 

As already shown by M08 and Belt et al. (2007) (at 
2 = 0), Eq. (|11|) does not provide a good fit to the data. 
In order to get a reasonable fit we need to add an explicit 
redshift dependence to the parameters a and /?, while in 
original the model proposed by AUgood et al. , the red- 
shift dependence of the shape parameter s was described by 
the quantity M*(2:). A second drawback of Eq. (|lll) can be 
clearly seen in figures [H [7] and [8] The convex-curved shape 
of the data deviates from a simple power-law behavior. In 
order to better model the data we modified Eq. (|ll|l into 



s{z, M) = Q(log(Mvir/[h-' Mq])* -f /?. 



(12) 



The values of the new fitting parameters a and /3 are re- 
ported in table 13] in the appendix. 

Finally we compare the inner (so. 3) to the outer shape 
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Figure 7. Mass and time evolution of the shape of dark matter 
haloes quantified via the p = 03/02 ratio. Points represent our 
data while the dashed line shows the fitting to equation 1121 to 
^ = 0. Similar results for the fitting are obtained for z > and 
for clarity are not shown in the plot. The color code is the same 
used in figure [T] and shows different rcdshifts. Values of the fit 
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Figure 8. Mass and time evolution of the shape of dark matter 
haloes quantified via the S0.3 = aa/ai ratio measured at 0.3-Rvir- 
Points represent our data while the dashed line shows the fitting 
to equation ll2l to 2: = 0. Similar results for the fitting are obtained 
for 2 > and for clarity are not shown in the plot. The color code 
is the same used in figure[T]and shows different redshifts. Values of 
the fit parameters for all redshifts can be found in the appendix. 



(s); we found that haloes are slightly aspherical in their cen- 
tral region (5-10% efTect) as already noticed by AUgood et 
al. (2006), who ascribed this difference to the higher num- 
ber of substructure in the central region with respect to the 
external one. However it should be kept in mind that the 
internal region is the most likely to be affected by baryonic 
physics (e.g. Kazantzidis et al. 2004) and hence the inclu- 
sion of dissipational effects is needed in order to drive more 
quantitative conclusions. 



5 HALO SPIN PARAMETER 

The relation between spin parameter (A') and halo mass 
is shown in figure [O] where for each redshift value we fit 
the data with a power law and in analogy with eq. [5] we 
denote the gradient and intercept in log space as a{z) and 
b respectively. The parameters of these fits are shown in 
the Appendix. For low redshift (z = 0; 0.23) we found our 
results to be consistent with the spin parameter being almost 
mass independent, in agreement with Maccio et al. (2007; 
2008) and Bett et al. 2007. At higher redshift {z > 0.5) a 
mass dependence of the spin starts to develop; with high 
mass haloes having on average a lower value of A'. This is 
consistent with the results of Knebe & Power (2008), who 
performed a detailed study of the spin-mass relation for z = 
1 and z — 10. 

The evolution of the slope of the spin-mass relation is 
more evident in figure [TUl where the values of the a{z) pa- 
rameter are shown (see eq:[5]). As can be seen in the figure, 
a evolves from being consistent with at z = to negative 
values for increasing redshift. Knebe & Power (2008) noted 
that this behavior is almost insensitive to the mass binning 



and details of the selection criteria for the relaxed halo pop- 
ulation, but they found that the number of particles per halo 
is an important parameter. In order to verify the stability 
of our results against the minimum number of particles per 
halo, we show in figure [10] the redshift evolution of the slope 
of the spin-mass relation for two different choices of iVvir, 
namely 500 and 1000 particles. As can be seen in the figure, 
the change in A^vir does change (slightly) the value of the 
fitted slope a{z), nevertheless this change is still within the 
statistical errors and does not affect the overall trend. 

It could be possible that the mass dependence we ob- 
serve could be partially due to the mixing of different box 
sizes (and hence haloes with different resolution). For this 
purpose we analyzed the spin-mass relation independently 
in different boxes (B90, B180 and B300). We found no signif- 
icant differences between our global analysis and the single 
box results. We can conclude that no numerical artifacts 
are affecting our results and that indeed spin parameter and 
halo mass are weakly correlated at high redshift. 

Finally, Fig. [TT] shows the distribution of spin parame- 
ters at redshifts 2 = 0, 1.12 and 2. At all redshifts the dis- 
tribution is well fitted by a log-normal distribution. 



P(A') = 



1 



A'V27rcr 



exp 



2cr2 



(13) 



with a — 0.57 and Aq — 0.031. This last quantity Aq shows 
a marginal mass dependence at high redshift as shown in 
figure |9] 
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Figure 9. Mass and time evolution of the spin parameter A'. 
(Top) Points are the median from our simulations while the solid 
lines represent the linear fitting to the data. Every line have been 
shifted by a constant factor of 0.1 from 2 = 2 and all of them show 
approximately the same mean value, nevertheless as z increases, 
a weak dependence of the spin on the mass of the haloes starts to 
become evident. The color code is the same used in figure [l] and 
shows different redshifts. (Bottom) Unshifted linear fitting to A'. 
The color code of each line is the same used in the plot on top. 
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Figure 10. Time evolution of the parameters of fitting a{z) and 
b{z) to equation |5] for the spin-mass relation for haloes with at 
least 500 (blue) and 1000 (red) particles inside the virial radius. 
Error bars show the error in the value of the parameter a given 
from the fitting. 
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and a = 0.57. 



6 DISCUSSION AND CONCLUSIONS 

In this work we present a detailed analysis of a large set 
of N-body simulations performed within a WMAP 5th year 
ACDM cosmology (Komatsu et al. 2009) . We study the re- 
lation between structural properties of dark matter haloes 
(concentration, spin and shape) with mass, and the evo- 
lution of such scaling relations with redshift. We span the 
entire mass range important for galaxy formation [10^° : 
W^^h~^ Mq] and a redshift range from z — to z = 2. 

We present results for "relaxed" haloes, defined accord- 
ing to the criteria suggested by Maccio et al. (2007). In our 
mass and redshift range the Cvir — Mvir relation always fol- 
lows a power law behavior. We confirmed that the redshift 
dependence of such relation is more complex than a sim- 
ple (1 + z)^^ scaling as proposed by Bullock et al. 2001a, 
with both the normalization and the slope of the relation 
changing with cosmic time. We also found that for increas- 
ing redshifts (z « 2) the power law behaviour seems to break 
in agreement with recent studies (e.g. Klypin et al. 2010). 
Thanks to our multiple box simulations we tested our results 
against resolution effects and find them to be stable once a 
sufficient large number of particles is used A'vir > 500. 

Recently two other works have addressed the topic of 
the evolution of the Cvir — Mvir relation, Zhao et al. (2009) 
and Klypin et al. (2010). When compared with the model 
proposed in Zhao et al. (2009) our results show a very good 
agreement at the low mass end (possibly due to the fact that 
both halo samples have more or less the same level of resolu- 
tion). For high masses we find a slightly higher difference but 
never exceeding few percent. The comparison with Klypin et 
al. is less straightforward since they used a different method 
to compute concentrations, based on the circular velocity of 
the halo instead of directly fitting the density profile. More- 
over they use all haloes in their simulation volume without 
any distinction between relaxed and unrelaxed. Our results 
are in qualitative agreement with the model proposed by 
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Klypin et al. for the evolution of the Cvir — Mvir relation but 
there are differences of the order of 8%. It is then interesting 
to ask ourselves if these differences arise from the different 
method used to compute Cvir. 

For this purpose we applied to our B90 box (which 
has roughly the same resolution as the Bolshoi simulation 
of Klypin et al. 2010) the method proposed by Klypin to 
compute the concentration, based on the relation between 
mass and maximum circular velocity. To be consistent with 
Klypin et al. 2010 we considered all haloes in our sample, 
without making any distinction between relaxed and unre- 
laxed. We find that the different methods to compute Cvir 
do not introduce any systematic bias, being perfectly con- 
sistent. Nevertheless we found that at fixed mass our haloes 
have a lower circular velocity compared to the Klypin et 
al. (2010) results. Taking into account the explicit depen- 
dence of the concentration on Vc we are keen to conclude 
that part of the difference in Cvir — Mvir relation could be due 
the slightly different values of the cosmological parameters 
and primordial power spectral index. 

Another interesting result regards the effect of including 
unrelaxed haloes. Klypin et al. (2010) did not make any 
attempts of removing unrelaxed haloes; this is in principle 
justified because their concentrations are obtained from an 
integral quantity (Vc) which is less affected by the dynamical 
status of the halo with respect to a differential one like p. 
On the other hand by applying Klypin's method to our B90 
box we found that the Vc — M relation has a systematic 
shift (« 6%) towards higher values when it is computed for 
"relaxed" haloes instead than for all haloes. This results in a 
systematic bias, towards lower concentrations, for the "all" 
sample with respect to the relaxed one. 

Finally Zhao et al. (2009) and Klypin et al. (2010) found 
that the concentration may have a minimum value close to 
3.5-4.5 at redshift close to 4. Unfortunately our simulations 
do not have enough resolution to get a statistical valid halo 
sample at such a high redshift, so we cannot confirm such 
a finding, even if we do see some evidence of a breaking of 
a simple power law behaviour for the Cvir — Mvir relation at 
z = 2. Let us stress one more time that the fitting formula we 
proposed in this work are valid only in the test redshift range 
[0-2] and should not be extrapolated at higher redshifts. 

In order to improve our understanding on the redshift 
evolution of the Cvir — Mvir relation we look at the individ- 
ual evolution with time of Ts and -Rvir. Both these length 
scales grow with decreasing redshift until a maximum is 
reached, then they start to decrease towards z = 0. There 
is a clear analogy between the collapse of a linear pertur- 
bation and the behaviour of rg and -Rvir. We found that we 
can model the evolution of the inner part of the halo (within 
Ts) as a decoupled spherical perturbation growing inside the 
central region of the halo. The temporal offset between the 
"turning points" of the perturbations associated with r^ and 
-Rvir is able to explain the observed redshift evolution of the 
Cvir — Mvir relation; which strongly deviates from a simple 
(1 + z)" scaling of the z = Q relation. As a final remark we 
would like to stress that our results refer to the dark mat- 
ter distrubution in the absence of a coUisional component. 
Although it has been shown that the inclusion of baryonic 
physics may affect the properties of the dark matter, it also 
known that the strength of this effect strongly depends on 
the implemented baryonic physics (Duffy et al. 2010, Gover- 



nato et al. 2010). The dark matter distribution in real haloes 
is still under debate and a comparison of pure dark matter 
results with observations should be preformed with extreme 
caution. 

We then investigate the mass and redshift dependence 
of the axis ratio (p and s) of dark matter haloes. Our re- 
sults are in agreement with previous studies and show that 
although on average haloes in our simulations are preferen- 
tially triaxial at all masses and redshifts, low mass haloes 
are more spherical than high mass ones and, at any mass, 
central regions are more aspherical than outer ones. We find 
a more complex evolution of the shape-mass relation with 
redshift with respect to the model proposed by AUgood et 
al. (2006). We propose a new fitting function (that deviates 
from a simple power law) that is able to reproduce our data 
in the whole redshift and mass range. 

Finally we studied the evolution of the spin parame- 
ter A'. At redshift zero we confirm the results of Maccio 
et al. (2007; 2008), who found the spin parameter to be 
mass independent. For increasing redshift there is evidence 
of a correlation between halo mass and spin: on average, 
more massive haloes have lower values of the spin parame- 
ter. This is in agreement with recent findings at very high 
redshift {z = 10) by Knebe & Power (2008). As already 
speculated by those authors, since disk sizes are to first or- 
der proportional to the spin parameter (Mo et al. 1998), a 
lower spin parameter for high mass haloes could make their 
central stellar and gaseous body more compact and hence 
allowing more efficient star formation. In addition this could 
affect the evolution of the size-mass relation of galaxy disks, 
which is usually modeled by assuming A is independent of 
mass and redshift (Mao, Mo & White 1998; Somerville et 
al. 2008; Firmani & Avila-Reese 2009; Button et al. 2010). 

However, in our simulations the differences in median 
spin parameters aX z = 2 compared to z; = are at most 
15%, and thus the consequences to observables such as the 
size-mass relation of disk galajcies are likely to be small. Fur- 
thermore, baryonic effects such as supernova driven outflows 
and inefficient cooling can modify spin parameters by much 
larger amounts (factors of ~ 2) (Button & van den Bosch 
2009). Thus it is unlikely that the mass dependence on halo 
spin that we observe in our simulations at 2; = 2 will have 
an unambiguous observational signature. 
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APPENDIX A: DETAILED VALUES FOR THE 
PARAMETERS OF THE DIFFERENT FITS 
PRESENTED IN THE MAIN BODY OF THE 
PAPER. 

We give in this appendix a complete set of tables with 
the value of all the fitting parameters presented in the pa- 
per. All fits to the data are performed for the functional 
form log('(/;) — alog(Mvir) + b where V can be either Cvir, 
A', s or p and where Mvir is in units of lO^'^/i^^M©. We 
stress on the validity of our results in the range of redshifts 
between and 2. Nevertheless individual fits are valid in 
the individual mass ranges [10^°/i"^ Mq, lO^'^h"^ M©] for 
z=0, [10^°h-'^MQ,~ 3.16"h-^M0] for z=0.23, 0.38 and 
0.56, [10^°/i"^ Mq, ~ 1.26"h-i Mq] for z=0.8, and 1.12 and 
[10^°/i"^ Mq, ~ 5.0^^h-^ Mq] for z=1.59, and 2.0. 

We present in table [2l just for completeness, the fitting 
of the shape paramaters to log(i/j) — alog(Mvir) + b for all 
masses and redshifts. In table[3]we present the fit to Eg. 1121 
to the shape parameters s,p and so. 3. 

This paper has been typeset from a T^j]X/ ffl^jX file prepared 
by the author. 
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Table 2. Values of the fit parameters for the data presented in 
the paper, log(cvir), log(A') as computed from haloes with at least 
500 and 1000 particles inside the virial radius and log(s) and 
log(p) as computed for haloes with more than 4000 particles inside 
the virial radius. All fits where done to the function log('(/') = 
alog(Mvir) + ^. N haloes is the total number of relaxed haloes 
used at that redshift to compute the mean and median values. 



Table 3. Values of the fit parameters for the data presented in 
the paper s, p and so. 3 as computed for haloes with more than 
1000 and 4000 particles inside the virial radius respectivelly. All 
fits where done to the function ip = Q[log(Mvir)]^ + /9. N haloes 
is the total number of relaxed haloes used at that redshift to 
compute the mean and median values. 



Redshift 


N haloef 


i a 


Aa 


b 


Ab 


Redshift 


N haloes 


a(xlO-'') 


Aa(xlO-'') 





A/3 




log( 


Cvir) (Nmin 


L = 500) 










s (Nmin = 


1000) 









23777 


-0.097 


0.002 


2.155 


0.021 


0.0 


12426 


-6.566 


2.760 


0.815 


0.008 


0.23 


22358 


-0.091 


0.002 


2.011 


0.027 


0.23 


11563 


-7.622 


2.818 


0.815 


0.008 


0.38 


20906 


-0.085 


0.001 


1.897 


0.018 


0.38 


10753 


-7.009 


1.185 


0.783 


0.003 


0.56 


19475 


-0.078 


0.001 


1.763 


0.013 


0.56 


10049 


-7.215 


1.484 


0.771 


0.004 


0.8 


17696 


-0.075 


0.002 


1.684 


0.030 


0.8 


8942 


-7.173 


3.320 


0.751 


0.008 


1.12 


14960 


-0.066 


0.001 


1.523 


0.017 


1.12 


7482 


-6.388 


3.945 


0.714 


0.010 


1.59 


11888 


-0.049 


0.003 


1.257 


0.032 


1.59 


5851 


-5.932 


3.760 


0.677 


0.008 


2.0 


9290 


-0.039 


0.005 


1.093 


0.055 


2.0 


4404 


-5.512 


3.936 


0.655 


0.009 




log 


(A') (Nmin 


= 500) 










p (Nmin = 


: 1000) 









23777 


-0.001 


0.004 


-1.497 


0.046 


0.0 


12426 


-6.345 


2.510 


0.955 


0.007 


0.23 


22358 


-0.000 


0.003 


-1.509 


0.038 


0.23 


11563 


-7.396 


7.872 


0.955 


0.023 


0.38 


20906 


-0.006 


0.004 


-1.429 


0.044 


0.38 


10753 


-6.452 


2.024 


0.923 


0.005 


0.56 
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-0.015 


0.005 


-1.308 
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0.56 


10049 


-6.802 


1.755 


0.916 
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0.8 
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-0.011 
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-1.369 


0.043 
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-7.539 


2.984 


0.916 
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1.12 
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-0.020 


0.004 
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1.12 
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-5.992 


5.600 
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0.014 


1.59 


11888 


-0.024 


0.004 


-1.220 


0.048 


1.59 


5851 


-6.299 


4.812 


0.845 


0.011 


2.0 


9290 


-0.030 
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0.064 


2.0 


4404 


-5.118 


7.149 


0.808 


0.016 




log( 


;A') (Nmin : 


= 1000) 










so. 3 (Nmin 


= 4000) 









12427 


-0.002 


0.005 


-1.485 


0.063 
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-6.464 


2.709 


0.760 


0.008 


0.23 
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-0.004 
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-1.452 
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0.23 
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3.404 


0.760 
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0.38 
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-0.010 
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-1.373 


0.066 


0.38 
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-6.776 


2.803 


0.722 


0.008 


0.56 
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-0.012 
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-1.337 
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0.56 
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-6.870 
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0.703 
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0.8 
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-0.013 
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-1.344 
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-6.444 


3.397 
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1.12 
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-0.016 
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-1.292 
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1.59 
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-0.030 


0.005 


-1.146 


0.067 


1.59 


1211 


-6.441 


7.297 


0.632 


0.020 


2.0 


4404 


-0.035 


0.011 


-1.091 


0.137 


2.0 


879 


-4.050 


0.138 


0.562 


0.035 




log 


(s) (Nmin = 


= 4000) 





















3102 


-0.039 


0.006 


0.386 


0.079 














0.23 


2785 


-0.033 


0.003 


0.294 


0.044 














0.38 


2523 


-0.036 


0.004 
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0.059 














0.56 


2281 


-0.035 


0.003 


0.303 


0.034 
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1972 


-0.038 


0.004 


0.333 


0.055 














1.12 


1588 


-0.033 


0.005 
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0.063 














1.59 


1211 
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0.005 
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0.059 














2.0 


879 
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0.010 


0.133 
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log 
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-0.051 


0.007 


0.443 


0.096 
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0.005 
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0.062 
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0.005 
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0.065 














0.56 


2281 
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0.004 
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0.051 














0.8 
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-0.044 
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0.298 


0.071 














1.12 


1588 


-0.039 


0.006 


0.224 


0.076 














1.59 


1211 


-0.034 


0.005 


0.145 


0.059 














2.0 


879 


-0.032 


0.012 


0.112 
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